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ABSTRACT 

We present the properties of the ensemble variability V for nearly 5,000 near-infrared active galactic 
nuclei (AGNs) selected from the catalog of Quasars and Active Galactic Nuclei (13th Edition) and 
the SDSS-DR7 quasar catalog. From three near-infrared point source catalogs; namely 2MASS, DE- 
NIS, and UKIDSS/LAS catalogs, we extract 2MASS-DENIS and 2MASS-UKIDSS counterparts for 
cataloged AGNs by cross-identification between catalogs. We further select variable AGNs based on 
an optimal criterion for selecting the variable sources. The sample objects are divided into subsets 
according to whether near-infrared light originates by optical emission or by near-infrared emission, in 
the rest-frame; and we examine the correlations of the ensemble variability with the rest-frame wave- 
length, redshift, luminosity, and rest-frame time lag. In addition, we also examine the correlations 
of variability amplitude with optical variability, radio intensity, and radio-to-optical flux ratio. The 
rest-frame optical variability of our samples shows negative correlations with luminosity and positive 
correlations with rest-frame time lag (i.e., the structure function, SF), and this result is consistent with 
previous analyses. However, no well-known negative correlation exists between the rest-frame wave- 
length and optical variability. This inconsistency might be due to a biased sampling of high-redshift 
AGNs. Near-infrared variability in the rest-frame is anticorrelated with the rest-frame wavelength, 
which is consistent with previous suggestions. However, correlations of near-infrared variability with 
luminosity and rest-frame time lag are the opposite of these correlations of the optical variability; 
that is, the near-infrared variability is positively correlated with luminosity but negatively correlated 
with the rest-frame time lag. Because these trends are qualitatively consistent with the properties 
of radio-loud quasars reported by some previous studies, most of our sample objects are probably 
radio-loud quasars. Finally, we also discuss the negative correlations seen in the near-infrared SFs. 
Subject headings: galaxies: active - quasars: general - catalogs 



1. INTRODUCTION 

Variability is a characteristic property of active galactic 
nuclei (AGNs). It is believed that the majority of AGNs 
exhibit rapid and apparently random variability over a 
wide wavelength band that ranges from x-ray to radio 
wavelengths. Characteristic variability timescales range 
from months to years. Variability is a powerful feature 
for constraining physical models and has the potential to 
probe the structure and dynamics of the central engine. 

It is widely accepted that accretion disk instabil- 
ities represent a promising model fo r explaining the 
varia bility properties of AGNs (e.g., iKawaguchi et al.1 
1998). In addition, several possible models have 
been proposed to explain the observed optical variabil- 
ity: m ultiple supernoy ae (iTerlevich et all 11992ft or star- 
burst (jTerlevich et al.lll992t lAretxaea et al.lll997f) close 
to the nucleus, star collisions jCourvoisier et al.l Q~996t 
iTorricelli-Ciamponi et all 20001) . and gravitationa l mi- 
crolensing (|Hawkinsl 11993 : iZackrisson et~aLl 120031) . A 
powerful approach for discerning these models is to ex- 
amine correlations between variability and a physical pa- 
rameter (such as time lag, luminosity, rest-frame wave- 
length, or redshift). These observational correlations 
have been associated with correlations predicted by phys- 
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ical models, and the appropriate model for the variability 
mechanism of AGNs has often been discussed. 

There are two main approaches to examine the vari- 
ability of AGNs. The first is to monitor a re latively small 
sample of AGNs over a period of time (e.g..|Garcia et al.l 
\19M INeugebauer fc Matthewsl[T999l : IWalsh et alJ l2009) . 
By using an observed light curve for an AGN, various 
variability properties for the object were studied. The 
second approach is to examine the ens emble variabil- 
ity V of a larger numb er of AGNs (e.g., Ide Vries et al.l 
2003: lBauer"etaII[2009l ). This approach can be accom- 
plished e ven if there are only two sampling epochs. For 
instance. iVanden Berk et al.1 (|2004j ) studied the ensemble 
variability of 25,000 quasars in the SDSS catalog using 
two sampling epochs and established several correlations 
between the optical variability and time lag, luminosity, 
rest-frame wavelength, a nd redshift. In the near-infrared, 
INeugebauer et al.l (|1989ft studied the near-infrared vari- 
ability of individual quasars using a sample of 108 opti- 
cally selected quasars. In add ition, near-infrared ensem- 
ble variability is reported by lEnva et "all (|2002ft . They 
presented correlations between near-infrared variability 
and several physical parameters by studying a sample of 
226 AGNs. Whereas numerous studies have been dedi- 
cated to examining the optical variability (especially to 
establish a correlation between optical variability and a 
physical parameter), a few studies have focused on the 
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near-infrared variability. 

In this paper, we present the optical and near-infrared 
ensemble variability in the rest-frame for a sample of 
about 5,000 AGNs, which are probably mostly quasars, 
by using archival three near-infrared catalogs and two 
AGN catalogs. The catalogs used in this paper are intro- 
duced in Section [5] In Section [31 we describe the method 
for extracting near-infrared counterparts for cataloged 
AGNs and the criterion for selecting near-infrared vari- 
able AGNs. In Section |4j we show correlations of the en- 
semble variability with rest-frame wavelength, redshift, 
luminosity, and rest-frame time lag. We also show the 
variability properties correlated with optical variability, 
radio intensity, and radio-to-optical flux ratio. Finally, 
we discuss both the intrinsic optical and near-infrared 
variability properties in Section [5] 

2. DATA 

2.1. 2MASS 

The Two Micron All Sky Survey (2MASS) 1 
(Skrutski e et al.l [200 61 is a project that observed 
99.998% of the whole sky at the J (1.25 ^m), H (1.65 
/xm), and Ks (2.16 /j,m) bands, from Mt. Hopkins, 
Arizona, USA (in the northern hemisphere) and from 
CTIO, Chile (in the southern hemisphere) between 
June, 1997 and February, 2001. Both the instruments 
are highly automated 1.3-m telescopes equipped with 
three-channel cameras, with each channel consisting 
of a 256 x 256 array of HgCdTe detectors. 2MASS 
obtained 4,121,439 FITS images (pixel size ~ 2"0) with 
an integration time of 7.8 s. The limiting magnitudes 
[signal-to-noise ratio S/N>10] are 15.8 (J), 15.1 (H), 
and 14.3 (K s ) mag at each band. The Point Source 
Catalog (PSC) was produced using these images and 
contains 470,992,970 sources. The images and the PSC 
are publicly available on the 2MASS web site. 

The 2MASS default magnitude in the PSC is measured 
by a point spread function (PSF) profile-fitting algo- 
rithm. If no source is detected, this is the 95% confidence 
upper limit derived from a 4" radius aperture measure- 
ment taken at the position of the source on the Atlas 
Image. The uncertainties in the profile-fit photometry 
take into account the PSF shape, using a measurement- 
noise model that includes the effects of read noise, Pois- 
son noise, and PSF error. 

2.2. DENIS 

The Deep Near Infrared Survey of the Southern sky 
(DENIS) 2 is a survey of the southern sky done between 
1995 and 2001 in two near-infrared bands (J: 1.25 /im and 
K s : 2.15 /im) and one optical band (I: 0.82 ^m). The 
instrument is a 1-m ground-based telescope at La Silla 
(Chile) with limiting magnitudes of 16.5, 14.0, and 18.5 
in the three bands, respectively. The positional accuracy 
is better than 1". The DENIS catalog consists of a set 
of ~ 350 million point sources in approximately 16,700 
deg 2 of the southern sky. The PSF and several kinds of 
aperture magnitudes are available in the DENIS catalog. 
The photometric quality is defined by quality flags, which 
are represented by a number between and 100. 

1 2MASS web site (http://www.ipac.caltech.edu/2mass/) 

2 http:/ /cdsweb.u-strasbg.fr/denis.html 



2.3. UKIDSS 

The United Kingdom Infrared Telescope (UKIRT) In- 
frared Deep Sky Survey (UKIDSS) 3 ([Lawrence et al.l 
I2007T) be gan in May 2005 and was planned to survey 
7,500 deg 2 of the northern sky, extending over both high 
and low galactic latitudes. The survey i nstrument is 
the UKIRT Wide Field Camera (WFCAM; ICasali et al.l 
120071 ). which has four 2048 x 2048 Rockwell devices, on 
the UKIRT in Hawaii. The pixel scale is 0"4. The 
pipeline processing and science archive are described in 
Irwin et al (2008) and Hambly et al (2008). 

The UKIDSS consists of five surveys: the Large Area 
Survey (LAS), Galactic Plane Survey (GPS), Galac- 
tic Clusters Survey (GCS), Deep Extragalactic Sur- 
vey (DXS), and Ultra Deep Survey (UDS). In this pa- 
per, we use the LAS data from UKIDSS DR6. The 
UKIDSS/LAS covers an area of -4,000 deg 2 at high 
galactic latitude in th e YJHK filter s . Th e photometric 
system is described in lHewett et al.l (j2006[ ). 

The UKIDSS/LAS point source catalog contains the 
aperture magnitude measured using several aperture 
sizes. The photometric errors are less than 0.02 mag 
in all filters, an d the global uniformit y is as good as the 
2MASS survey (jLawrence et al.ll2007l ). 

2.4. Sample AGNs 

Sample AGNs are extracted from two catalogs: the 
cat alog of Quasars and Active Galactic Nuclei (13th 
Ed; iVeron-Cettv fc Veronll2010t hereafter OA) an d the 
SDSS-DR7 quasar catalog ([Schneider et al.l 120101 here- 
after SQ). The QA catalog includes 133,336 quasars, 
1,374 BL Lac objects, and 34,231 active galaxies (includ- 
ing 15,627 Seyfert 1) with positional accuracy superior to 
1"0. The SQ catalog includes 105,783 spectroscopically 
confirmed quasars in an area covering ~9,380 deg 2 , and 
has a positional accuracy better than 0"1 rms per coor- 
dinate. 

TABLE 1 

Number of cataloged AGNs having near-infrared 
counterparts. 





DENIS 


UKIDSS 




J K 


J K 


QA 


2750 1159 


3832 3540 


SQ 


551 110 


2484 1945 



3. SELECTION OF NEAR-INFRARED VARIABLE 

AGNS 

3.1. Extraction of Near-infrared Counterparts for 
Cataloged AGNs 

We first extracted a 2MASS counterpart for each cata- 
loged AGN by cross-identifying the 2MASS PSC with the 
QA and SQ catalogs. Each 2MASS counterpart was fur- 
ther cross-identified with the DENIS and UKIDSS/LAS 
catalogs, and we extracted near-infrared magnitudes for 
different epochs. We only use sources having sufficient 
photometric quality in either the J or K bands: 2MASS 
photometric quality flags superior to C [corresponding to 

3 http://www.ukidss.org/ 
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Fig. 1. — Magnitude versus magnitude diflercncc between 2MASS and either the DENIS or UKIDSS catalogs, and histograms with 
respect to magnitude difference. In each panel, the upper figure is for normal stars (black dots) and the bottom figure is for extracted 
AGNs (red dots). The standard deviation for magnitude differences is shown in the upper right of each histogram. 



S/N > 5 and a corrected band photometry uncertainty 
(cmsig)< 0.21714], DENIS quality flags superior to 90, 
and UKIDSS/LAS photometry having pNoise (proba- 
bility that the source is noise)< 0.1, ErrBits= 0, and 
ppErrBits= 0. The ErrBits (ppErrBits) indicates pro- 
cessing (post-processing) error quality bit flags, with zero 
being the best quality. Note that all magnitudes at K s 
(2MASS or DENIS) and K (UKIDSS) bands are denoted 
by K, throughout this paper. The number of the cata- 
loged AGNs having near-infrared counterparts is sum- 
marized in Table [T] 

3.2. Zero Point of Magnitude Difference 

To examine the near-infrared variability of the cata- 
loged AGNs, we use the PSF magnitudes in both the 
2MASS and DENIS catalogs. Because there are no PSF 
photometry measurements in the UKIDSS/LAS catalog, 
we use average aperture magnitudes between the 2. "8 
and 5. "7 apertures. Because we focus on the variability 
of AGNs by comparing the magnitudes listed in differ- 
ent catalogs, we should take into account the differences 
in the photometric systems between different catalogs. 
To examine a photometric difference, we first investi- 
gate magnitude differences of normal stars. We define 
all the magnitude differences to be DENIS or UKIDSS 
magnitudes subtracted from 2MASS magnitudes. We ex- 
tracted a sample of 2MASS stars based on near-in frared 
colors of normal stars from IBessell fc Brett] JT988). The 
DENIS and UKIDSS magnitudes for these stars were 
derived by catalog cross-identification. We draw his- 
tograms for the magnitude differences of normal stars, 
which show Gaussian distributions with small but non- 
zero peak values. However, AGNs with near-infrared 
counterparts have distributions similar to those of nor- 



mal stars but have larger peak values than normal stars. 
This might be because the AGNs are point sources but 
are barely more diffuse than the normal stars. Therefore, 
we corrected the magnitude differences so that the peak 
AGN values are coincident with 0. Here, we redefine the 
amplitude of variability for the i-th object as follows: 

A?7lj = 77li,2MASS — m i, DENIS or UKIDSS ~\~ ^ TO cor j (l) 

where Am„ is the correction term for the magnitude 
difference. We calculated correction terms at both the J 
and K bands for both the DENIS and UKIDSS catalogs 
(i.e., we derived four correction terms) and found the 
values to be -0.18 (DENIS- J), -0.12 (UKIDSS- J), -0.14 
(DENIS-K), and -0.10 (UKIDSS-K). 

Figure [T] shows magnitude versus corrected magni- 
tude difference, and histograms for each case. We also 
corrected the magnitude differences for normal stars, 
with the correction terms of -0.02 (DENIS- J), -0.02 
(UKIDSS- J), -0.03 (DENIS-K), and 0.0 (UKIDSS-K). 
As mentioned above, no significant difference in the his- 
togram shape between normal stars and extracted AGNs 
can be seen. The near-infrared counterparts for the ex- 
tracted AGNs are point sources because we extracted the 
counterparts from point source catalogs. Consequently, 
the similarity in the shapes is a natural result. Fur- 
thermore, the fact that the extracted AGNs are point 
sources sugg ests that mo s t of t he AGNs are expected to 
be quasars. ISesar et al.l (|2007l ) estimated that at least 
90% of quasars are variable at the 0.03 mag level. In all 
cases, the standard deviations of magnitude difference for 
the AGNs are larger than those for normal stars. This 
implies that most of the AGNs have a more-or-less vari- 
able nature at near-infrared wavelengths. 




Fig. 2. — Magnitude versus magnitude difference for the variable (red solid circles) and non-variable (blue crosses) AGNs extracted with 
our selection criterion. The curves show la confidence interval for magnitude differences of the normal stars, which are derived using the 
sample stars in Figure [T] 
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Fig. 3. — Redshift versus absolute magnitude. The left panel is the sample in the J-band (blue), and the right panel is the sample in the 
K-band (red). 



TABLE 2 

Number of the selected near-infrared variable AGNs. A 

SUBSET OF VARIABLE AGNS IS LISTED IN BOTH THE QA AND SQ 
CATALOGS. 





DENIS 


UKIDSS 




J K 


J K 


QA 


554 426 


1731 1497 


SQ 


105 55 


1524 1073 


Total 


592 445 


2316 1920 


3.3 


Sample 


Selection 



The three near-infrared catalogs contain photometric 
errors. Taking into account these photometric errors, 
we selected optimal variable AGNs. For this paper, we 
regard AGNs having magnitude differences larger than 
the sum of photometric errors as near-infrared variable 
AGNs. In other words, we established the following cri- 
terion for selecting variable AGNs: 

A771j > (Xj^MASS "I" other 7 (2) 

where c^mass is the 2MASS photometric error of the 
i-th object and <7i. othcr is either the DENIS or UKIDSS 
photometric error of the i-th object. The selection of 
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Fig. 4. — Variability as a function of rest-frame wavelength. Colors indicate either J-selection (blue) or K-selection (red). Solid squares 
and circles show samples from DENIS and samples from UKIDSS, respectively. Three different symbol sizes represent the number of data 
including in each data point. The largest size contains more than 9 data, the middle contains more than 4 and less than 10 data, and the 
smallest contains less than 5 data. 



variable AGNs is independently performed at both the 
J and K bands. The number of near-infrared variable 
AGNs is summarized in Tabled Figure [5] shows the vari- 
able and non- variable AGNs on magnitude versus mag- 
nitude difference diagrams, with curves showing la con- 
fidence interval for magnitude difference of the normal 
stars in Figure [T] As seen in the figure, very few sources 
with magnitude differences smaller than la level are se- 
lected as variable AGNs. Throughout this paper, we call 
the sample of variable AGNs selected in the J band (K 
band) as the J selection (K selection) or J-selected (PC- 
selected) sample and call the sample derived from the 
DENIS (UKIDSS) catalog in the J selection (K selec- 
tion) the DENIS (UKIDSS)-J (K) sample. 

Using the redshifts listed in the AGN catalogs, we 
calculated the absolute magnitudes for variable AGNs. 
Assumed cosmological constants are Hq = 71 km s _1 
Mpc-\ Q M = 0.3, and fl A = 0.7. We applied K- 
corrections based on a power-law spectrum with a = 0.5. 
Redshift versus absolute magnitude relations are shown 
in Figure [3] 



4. ANALYSIS 

4.1. Ensemble Variability 

Since we have only two sampling epochs for each AGN, 
the magnitude difference does not directly reflect the 
characteristic amplitude of intrinsic variability. In this 
case, the ensemble variability can be useful to examine 
the variability of AGN ensembles. 



Here, we define the ensemble variability as 



N 9 



V = 



N 



where N is the number of variable AGNs, and 

2 _ 2 , 
°i "'i^MASS ' ® 

The error in V is represented by 



2 

Z.othc 



(3) 
(4) 



a v 



1 

2V 



\ 



(Ef 



Am 



EN 2 
i °l 



2a? 



N 3 



N 2 



( 5 ) 

where try is the measurement uncertainty. These are 
essent ially the same definitions as used by lEnva et"aLl 

(pool . 

We derive the ensemble variability for each physical 
parameter bin (i.e., rest-frame wavelength, redshift, ab- 
solute magnitude, and rest- frame time lag). Accordingly, 
the ensemble variability will be a function of a physical 
value. Below, we examine the correlations between the 
ensemble variability and the physical parameters. 

4.2. Rest-frame Wavelength 

Quasar variability is known to be related to the rest- 
frame wavelength. In the optical, many previous stud- 
ies have demonstrated that the varia bility decreases 
with i ncreasing rest-frame wavelength ([ Giallong o et al 
119911: iCristiani etldl 119961: iCid Fernandes et all |19& " 
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Fig. 5. — Variability as a function of redshift. The symbols are the same as in Figure|4] 



ITrevese fe Vagnettil 12001 IVanden Berk et al.ll200l . In 
addition, anticorrelation has been establishe d to span at 



least f rom the ultraviolet to near-i n frared (|Cutri et al 
19851 iPaltani fc Couryoisierl 119941: Idi Clemente et al 
199a iCristiani et all Il997t ITrevese fc Vagnettil l2002t 



Meusinger et al.ll201lh . 

Figure Q] shows variability as a function of rest-frame 
wavelength. For all the samples, a clear trend in the 
ensemble variability is to decrease with increasing rest- 
frame wavelength in the near-infrared (i.e., A roat > 0.8). 
Although the data points for the DENIS-K sample are 
unreliable in the range of A roat < 1.5 because of low statis- 
tics, an anticorrelation over at least A rost ~ 1.5 is clearly 
evident. This general trend suggests that the anticorre- 
lation between variability and rest-wavelength is estab- 
lished not only in the optical but also in the near-infrared. 

Conversely, in the optical (i.e., A rcst < 0.8), variabil- 
ity appears to increase with increasing rest-frame wave- 
length. A turnover can be seen around A rGBt = 0.8 in the 
DENIS-J sample. This is significant because the Spear- 
man rank-correlation coefficient is 0.891 (—0.576) for the 
positive (negative) slope in A rcst < 0.8 (A rcst > 0.8), with 
the significant at 99% (91%) confidence level. Only the 
UKIDSS-J sample shows a moderately flat curve. These 
trends in the optical are inconsistent with previous stud- 
ies. A possible cause of these conflicting results is bi- 
ased sampling of high-rcdshift AGNs. AGNs with higher 
redshifts become more difficult to detect because their 
apparent magnitudes should be more faint. Even if an 
AGN could be detected at its brighter phase, the mag- 
nitude is expected to be near the limiting magnitudes 
of the near-infrared catalogs. In most cases, the faint- 



phase magnitude of the AGN is expected to be fainter 
than the limiting magnitudes. Because we selected vari- 
able AGNs by comparing only double-epoch magnitudes, 
we will not be able to, in such a case, extract the AGN 
as a variable source. In particular, AGNs having larger 
variability tend to fall into this situation and therefore 
cannot be selected as variable AGNs. Accordingly, we 
could not extract variable AGNs having large variability 
amplitudes even if such AGNs exist, and most variable 
AGNs we could extract should have relatively small vari- 
ability amplitudes. Because AGNs with higher redshifts 
are detected at shorter wavelengths in the rest-frame for 
a fixed passband in the observer-frame, this biased sam- 
pling would especially affect the correlation at shorter 
optical wavelengths. 

Indeed, the steeper slope for the DENIS-J sample than 
that for the UKIDSS-J sample can be interpreted as 
the difference between limiting magnitudes. The DE- 
NIS cannot detect fainter AGNs in the same way as the 
UKIDSS, because the limiting magnitudes of the DE- 
NIS are brighter than that of the UKIDSS. As a result, 
at shorter rest-frame wavelengths, the number of high- 
rcdshift AGNs showing large variability in the DENIS-J 
sample becomes smaller than that in the UKIDSS-J sam- 
ple due to the above-mentioned biased sampling. This 
might cause the steeper slope in the DENIS-J sample. 

4.3. Redshift 

Correlations between redshift and variability are shown 
in Figure [SJ In the UKIDSS-J sample, variability in- 
creases with increasing redshift below z ~ 0.8. The 
DENIS-J sample also shows an increase similar to that 
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Fig. 6. — Variability as a function of near-infrared absolute magnitude. The symbols are the same as for Figure [4] The J-selected 
(K-selected) samples are separated at z = 0.5 (z = 1.6). 



of the UKIDSS-J sample. However, the slopes appear to 
change significantly around z — 0.6-0.8. Above z ~ 0.8, 
there seems to be an anticorrelation in both samples, al- 
though many data points are unreliable because of the 
small number of objects contributing to the data (this 
is especially true for the DENIS- J sample). The J- band 
wavelength in the observer-frame corresponds to a wave- 
length of 0.70-0.78 /im in the rest-frame at z = 0.6-0.8. 
This is a near boundary wavelength between the optical 
and the near-infrared regions. 

The K-selected samples also show a property similar 
to the J-selected samples. The UKIDSS-K sample shows 
a positive correlation below z ~ 2. In the rest- frame 
at this redshift, the K-band wavelength in the observer- 
frame also corresponds to a near boundary wavelength 
between the optical and the near-infrared. Beyond z ~ 2, 
there seems to be a negative correlation. No trend can be 
found in the DENIS-K sample because most data points 
are scattered because of the statistically small number of 
d ata at each p oint. 

lEnva et "all (|2002l ) reported that the sample of radio- 
loud AGNs shows a positive correlation between near- 
infrared variability and redshift, whereas the sample of 
radio-quiet AGNs shows a negative correlation. Accord- 
ing to their correlations, our sample AGNs exhibit a 
property similar to their radio-loud AGNs. 

Note that the correlations of variability with 
redshift can be interpreted as an anticorrelation 
between wavelength and variability (|Cristiani et all 
1991 iCid Fernandes et aLlll996tldi~Clemente et allll996l: 
Cristiani et al.l 119971 ). In other words, for a fixed pass- 
band in the observer-frame, AGNs with larger redshifts 
are detected at shorter wavelengths, which have system- 
atically large fluctuations. In this case, the correlation 
of variability with redshift does not reflect evolution- 
ary effects. Therefore, the variations in the slope near 
the boundary wavelength between the optical and near- 
infrared should be intrinsic to the variations in the slope 



mentioned in Section l4~2l 

4.4. Absolute Magnitude 

As shown above, the intrinsic variability properties in 
the optical appear to be different from those in the near- 
infrared. Consequently, in what follows, we divide our J- 
selected (K-selected) samples into two subsets according 
to whether an object's redshift is z < 0.5 or z > 0.5 (z < 
1.6 or z > 1.6). These boundary redshifts correspond 
to a rest-frame wavelength of about 0.8 /an, which is a 
near boundary wavelength between the optical and near- 
infrared. When the redshift of an object is larger than 
the boundary redshift (i.e., z = 0.5 or z = 1.6 for J- or 
K-selected samples, respectively), the variability of the 
object represents the intrinsic optical variability. 

Figure [6] shows correlations between absolute magni- 
tude (luminosity) and ensemble variability. For "all" of 
the J-selected samples (see bottom-row panels in Fig- 
ure [6]), a broad peak appears roughly at M 3 ~ —27. In 
other words, the ensemble variability increases with de- 
creasing absolute magnitude (i.e., negative correlation) 
for Mj > —27, whereas there is positive correlation for 
Mj < —27. As can be seen from the middle- and top- 
row panels in Figure [5J the slope variations are caused 
by a difference in variability between the optical and the 
near-infrared. The K-selected samples also show simi- 
lar trends; that is, a negative correlation appears in the 
near-infrared (i.e., < z < 1.6) samples whereas a pos- 
itive correlation appears in the optical (i.e., z > 1.6) 
samples. 

Various studies have confirmed that the optical vari- 
ability of quasars is anticorrelated with luminosity 
(lUomoto et all 119761 : iHook et all 119941: iCristiani et all 
119961: IGiveon et all I1999L IVanden Berk et all I2004D . 
which indicates that more-luminous quasars exhibit 
smaller variability. The optical variability of our sam- 
ples is consistent with these analyses. 

In Poissonian models, the following relationship is ex- 
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Fig. 7. — Structure functions for each sample. The symbols are same as Figure [4] The J-selected (K-selected) samples are separated at 
z = 0.5 (z = 1.6). 



pcctcd between fractional optical variability AL/L and 
mean luminosity L: 

— (X (6) 

where P = \ (|Cid Fernandes et al.ll2000h . Equation @ 
is equivalent to 

log(V) = + K, (7) 

where M is the absolute magnitude and K is a con- 
stant. Using this relationship, we measure, by weighted 
least square fits, /3 = 0.096 (DENIS- J), 0.080'(UKIDSS- 
J), 0.396 (DENIS-K), and 0.167 (UKIDSS-K). Only the 
DENIS-K sample has a steeper slope, which is consider- 
ably steeper than the other three samples. Despite this 
result, we expect the DENIS-K slope to be less reliable 
than the others because it is calculated from only two 
data points. The slopes of the other three samples are 
less than 0.2, which is inconsistent with the Poissonian 
prediction. Additionally, the three slopes are also less 
than such as j3 = 0.246 ±0.005 (iVanden Berk et al.l l200l 
and p = 0.205 ± 0.002 (jBauer et al.ll2009ft . We consider 
it highly probable that the more-level slopes are due to 
the same cause as mentioned in Section 14.21 namely, to 
a biased sampling of large-redshift AGNs. 

In contrast with the optical correlations, the near- 
infrared variability is negatively correlated with abso- 
lute magnitude for the whole sample. The slope of 
each sample is -0.062 (DENIS- J), -0.156 (UKIDSS-J), 
-0.093 (DENIS-K), and -0.120 (UKIDSS-K). In addi- 
tion, these negative correlations are also qualitatively 
co nsistent with the t rend of radio-loud AGNs reported 
bv lEnvaet all ([2002T ). 

The near-infr ared correlations can be e xplained by the 
calculations of iKawaguchi fc Moril (|201lD . which predict 
that thinner tori have larger near-infrared variability am- 
plitudes (see Figure 10 in their paper). Assuming more 



luminous A GNs have thinner tori [i.e., the receding torus 
model (e.g., lLawrencdfr991|) ]. their model implies that 
more-luminous AGNs have larger near-infrared variabil- 
ity. 

4.5. Structure Function 

The ensemble variability as a function of rest-frame 
time lag is identical to the structure function (SF). The 
SF has been used to study both individual quasars and 
quasar ensembles and is useful for understanding the 
timescales of the relevant variability processes. 

The logarithmic slope (7) of the SF can be an indicator 
for discriminating physical models and is described as 
follows: 

Uocr 7 , (8) 

where r is the rest-frame time lag with the time di- 
lation effect caused by cosmic expansion, taken into 
account. The rest- frame time lag is given by r = 
At/ (1+z), where At is t he time lag in the observer-frame. 
IKawaguchi et al.1 (jl998) derived 7 values by constructing 
theoretical models. Their results are 7 ~ 0.74-0.90 in the 
starburst model and 7 ~ 0.41-0.49 in the disk instabil- 
ity mode l . Add itionally, by using a microlensing model, 
Hawkins! (|2002l ) predicted 7 = 0.25 ± 0.03. 

Figure shows the SF for each sample. As done in 
Section 14. A\ the samples are divided into subsets ac- 
cording to the redshift. In the optical, the ensemble 
variability tends to increase with increasing rest-frame 
time lag for the UKIDSS-J and UKIDSS-K samples. 
The logarithmic slopes of these samples are computed 
by weighted least square fits to be 7 = 0.13 (UKIDSS- 
J) and 7 = 0.66 (UKIDSS-K), respectively. Several 
studies have shown that the SF in the optical increases 
steadily for r < 2-3 y r and becomes relatively con- 
stant for larger time lags dHook et al.ll994t iTrevese et al.l 
11994 iCristiani et"al1 [1990). Accordingly, we attribute 
the more-level slope in the UKIDSS-J sample to the dif- 
ference in time ranges. In other words, the slope of 
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Fi g. 8. — Optical variabil ity (SDSS r-band) versus near-infrared variability amplitude. The optical variability for each object is taken 
from Mcusingcr ct al. (2011). 



the SF of the UKIDSS-J sample is expected to be less 
than that of the UKIDSS-K sample because a majority 
of the data in the UKIDSS-J sample have time lags longer 
than t = 2-3 yr whereas a majority of the data in the 
UKIDSS-K sample have time lags shorter than r = 2-3 
yr. Meanwhile, the results for the DENIS samples are 
inconsistent with those from both the UKIDSS samples 
and previous studies. Although the ensemble variability 
of the DENIS-J sample decreases with increasing rest- 
frame time lag, the number of objects in the DENIS-J 
sample is considerably smaller than that in the UKIDSS- 
J sample. The DENIS-K sample, with z > 1.6, has a no- 
tably smaller number of objects than the other samples. 
Therefore, we expect the results derived from the DENIS 
samples to be less reliable than those derived from the 
UKIDSS samples. 

In contrast, in the near-infrared, the ensemble vari- 
ability of all samples appears to be anticorrelated with 
the rest-frame time lag. The logarithmic slopes 7 are 
-0.002 (DENIS-J), -0.82 (UKIDSS-J), -0.14 (DENIS- 
K), and -0.67 (UKIDSS-K). Even the DENIS samples 
show negative correlations with the rest-frame time lag, 
despite their being relatively unreliable because of the 
same reason mentioned in the above. If these trends are 
real, then Figure [7] shows that the decrease is moderate 
at shorter time lags (r < 1000) whereas the decrease is 
steep at larger time lags (r > 1000). We further discuss 
the near- infrared SFs in Section [5j 



4.6. Multiwavelength Properties 
4.6.1. Optical Variability 

IMeusinger et al.l (|2011l ) studied the UV/optical ensem- 
ble properties of over 9,000 quasars from the SDSS Stripe 
82, and confirmed the well-known dependence of optical 
variability on time lag, luminosity, wavelength, and red- 
shift. The complete table of their quasar sample is avail- 
able at the CDS, which contains optical variability repre- 
sented by the variability indicator in their paper. Using 
their variability indicators, we examine here a correlation 
between optical variability and near-infrared variability. 

We first cross-identified our sample objects with their 
sample quasars and then extracted optical variability for 
each cross-identified object. We cross-identified 25 ob- 
jects for DENIS-J, 238 for UKIDSS-J, 22 for DENIS-K, 
and 173 for UKIDSS-K. Figure |S] shows the optical vari- 
ability at the SDSS r-band versus near-infrared variabil- 
ity amplitude. There seem to be positive correlations in 
the J-selected samples but negative correlations in the K- 
selected samples. The Spearman rank-correlation coeffi- 
cients are 0.113 (DENIS-J), 0.058 (UKIDSS-J), -0.107 
(DENIS-K), and -0.076 (UKIDSS-K) with the signifi- 
cant at 41%, 62%, 36%, and 68% confidence levels, re- 
spectively. Although these levels are insufficient, there 
might be a weak correlation in each sample. Note that 
the correlations for the J-selected samples are opposite 
those for the K-selected samples. 
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Fig. 9. — Radio flux versus near-infrared variability amplitude. Radio flux values taken from either the NVSS or FIRST survey. 



Because the near-infrared variability of our samples is 
simply the magnitude difference between double-epoch 
magnitudes, it does not directly represent the character- 
istic amplitude of the intrinsic variability. In addition, 
there are a statistically small number of objects in the 
samples. Consequently, the correlations must be con- 
firmed by using a sample with a larger number of objects. 

4.6.2. Radio 

Several studies have reported that radio-loud quasars 
tend to have larger optical variabi l ity amplitude than 
radio -quiet quasars (|Pica fc Sm ith 19831; iGarcia et al.l 
1999). However, some recent studies have shown that 
there is no significant dependence of optical variabil- 
ity amplitud e on ra dio loudness dHelfand et all 120011 : 
iBauer et all 120091 ) . iVanden Berk et al.l (|2004T )report 
that radio-loud quasars are more optically variable than 
radio-quiet quasars whereas there is no significant dif- 
ference in optical variability between radio-detected and 
ra dio-undetected q u asars . These trends were confirmed 
by iRengstorf et all (|2006f ) who used 933 known quasars 
in the QUEST Variability Survey. 

The SQ and QA catalogs contain radio flux extracted 
from the Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST) survey and either the National Ra- 
dio Astronomy Observatory Very Large Array Sky Sur- 
vey (NVSS) or FIRST surveys, respectively. DENIS-J 
has 209 objects with radio flux data, UKIDSS-J has 404, 
DENIS-K has 179, and UKIDSS-K has 353. 



We now review the relations between radio intensity 
and near-infrared variability. Figure [5] shows radio flux 
at 20 cm versus near-infrared variability amplitude. The 
radio flux is represented by the unit Jy. Because the 
SQ catalog provides 20 cm radio flux in the form of AB 
magnitude, we converted the radio magnitudes into radio 
flux in the unit of Jy. Each sample shows a weak cor- 
relation; that is, the near-infrared variability gradually 
increases with increasing radio intensity. The Spearman 
rank-correlation coefficients for the samples are 0.346 
(DENIS-J), 0.242 (UKIDSS-J), 0.235 (DENIS-K), and 
0.270 (UKIDSS-K). Each of these values has a signifi- 
cance level below 1% under the null hypothesis that vari- 
ability and radio intensity are not correlated. Note that 
we ruled out the sources having F 20 cm > 10, because 
they appear to be saturated as seen in the figure for the 
UKIDSS-K sample. 

We notice that positive correlations appear to be 
stronger above 0.01-0.1 Jy. Several studies have re- 
ported that radio-loud quasars tend to be more vari- 
able than radio-quiet quasars, although the optical vari- 
ability of simple radio-detected quasars is not signifi- 
cantly different from tha t of radio-undetected quasars 
(jVanden Berk et al.l [200l IRengstorf etail [2006]) . Our 
samples confirm these trends in the correlation between 
radio intensity and near-infrared variability amplitude. 

Moreover, there appears to be a positive correlation be- 
tween radio-loudness and near-infrared variability. Fig- 
ure [TU] shows the radio-to-optical flux ratio versus near- 
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Fig. 10. — Radio-to-optical flux ratio (F 2 o cm/Fu) versus near-infrared variability amplitude for variable AGNs selected from the SQ 
catalog. Both radio flux at 20 cm and optical flux at the SDSS u-band are taken from the SQ catalog. The left (right) panel shows the 
J-selected (K-selected) sample selected from the SQ catalog. The solid squares and open circles represent the DENIS and UKIDSS samples, 
respectively. 



infrared variability amplitude. We define here the radio- 
to-optical ratio as the ratio of radio flux at 20 cm to 
the optical flux at the SDSS u-band. Both radio and 
optical fluxes are taken from the SQ catalog. We find 
positive correlations in Figure 1101 which are significant 
because the Spearman rank-correlation coefficients are 
0.407 (DENIS- J), 0.230 (UKIDSS- J), 0.570 (DENIS-K), 
and 0.134 (UKIDSS-K) with the significant at 96%, 99%, 
97%, and 92% confidence levels, respectively. 

Theoretical calculations for near-infrared variability 
of AGNs predict larger variability a mplitude for orien- 
tatio ns closer to the pole-on view (Kawagu chi fc Mori 
2011). If the radio emission of our AGN sample origi- 
nates in (short or small) jets, then the radio intensity is 
likely to be correlated with the viewing angle in the sense 
that intense emission is associated to nearly pole-on view. 
Thus, the correlation between the near-infrared variabil- 
ity amplitude and the radio intensity that we discovered 
in this study is consistent with this prediction. 

5. DISCUSSION 
5.1. Optical Variability 

The intrinsic optical variability in both the DENIS 
and UKIDSS samples is inconsistent with the well-known 
negative correlation between rest-frame wavelength and 
variability. As mentioned in Section 14.21 this is prob- 
ably due to a biased sampling of high-redshift AGNs. 
However, the optical variability of the UKIDSS samples 
exhibit positive correlation with both the absolute mag- 
nitude and rest-frame time lag, which is consistent with 
previous analyses. On the other hand, although the ab- 
solute magnitude of the DENIS samples exhibit positive 
correlation with optical variability, the SF appears to 
be inconsistent with previous analyses. Therefore, com- 
pared with the DENIS samples, the UKIDSS samples 
strongly reflect the intrinsic properties of the optical vari- 
ability for AGNs. 

We attribute the result that the variability properties 
of the UKIDSS samples are more consistent than those 
of the DENIS samples to the following reasons: First, 
the UKIDSS samples contain significantly more number 
of objects (over three times more; see Table [2]) than the 
DENIS samples contain, and a statistically small number 



of objects increases errors. Second, photometric errors of 
the UKIDSS samples are roughly five times lesser than 
those of the DENIS samples. Consequently, both statis- 
tically and photometrically, the UKIDSS samples should 
be more reliable than the DENIS samples. 

5.2. Near- infrared Variability 

The near-infrared variability is anticorrelated with 
rest-frame wavelength, which is consistent with previous 
suggestions. Meanwhile, the near-infrared variability is 
anticorrelated with both absolute magnitude and rest- 
frame time lag, which are the opposite of these correla- 
tions of the optical variability. As mentioned in Section 
14.41 a subset of these correlatio ns are consis t ent w ith the 
radio-loud AGNs reported by lEnva et"all ((2002). De- 
spite the different optical properties between the DENIS 
and UKIDSS samples, both show similar near-infrared 
variability properties. 

We focus here on the near-infrared SFs of our sam- 
ples. The difference between the optical SF and the 
near-infrared SF is evident from Figure [7] Contrary to 
the optical SFs, the near-infrared SFs exhibit negative 
correlations. In particular, the UKIDSS samples show 
strong anticorrelations over ^1000 days. Note that the 
first-order SF is related to the a utocorrelation function 
as follows (jSimonetti et al.lfl985f ): 



V(t) = 2a 2 [l-p(r)}, 



(9) 



where a is the variance of a process and p(r) is the au- 
tocorrelation function with p(0) = 1. In this scheme, 
a negative SF slope indicates that the autocorrelation 
function increases with increasing rest-frame time lag. 
Accordingly, our samples imply that the near-infrared 
variability exhibits more periodic fluctuations on a longer 
time scale (up to ~10 yr). 

This anticorrelation can be seen in the SFs 
of several quasars in INeugebauer et al.l (|1989D and 
INeueebauer fe Matthews! (|1999f ). Several quasars in their 
samples show SFs that decrease over a period of 3-5 yr. 
In particular, the SF for quasar PG 1226+023 is appar- 
ent. O ur UKIDSS samples are c o nsiste nt with this prop- 
erty. Ncugeb auer fc Matthewsl (|1999l ) also insist that 
there is "periodic" (sinusoidal-like) behavior over nearly 
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a decade in the J, H, and K bands and suggest that this 
periodic behavior is a feature of the nonthermal engine 
of the radio-loud quasars. The peak near 3-5 yr that is 
seen in quasar PG 1226+023 appears to be found in the 
UKIDSS-K sample. Therefore, the SFs for the UKIDSS 
samples also show a periodic property similar to that 
found for radio-loud quasars, which is consistent with 
the analyses in Sections 14.31 and 14.41 

6. CONCLUSION 

We have examined the ensemble variability of over 
5,000 near-infrared luminous AGNs extracted from two 
AGN catalogs. The sample AGNs are divided into sub- 
sets according to whether, in the rest-frame, the near- 
infrared emission or the optical emission is at the origin of 
the near-infrared light, and we have demonstrated the de- 
pendence of intrinsic optical/near-infrared variability on 
rest-frame wavelength, redshift, absolute magnitude, and 
rest- frame time lag (i.e., structure function). The corre- 
lations of the optical variability with absolute magnitude 
and rest-frame time lag are qualitatively consistent with 
previous analyses, whereas the optical variability is in- 
consistent with the well-known negative correlation of 
optical variability with rest-frame wavelength. The in- 
consistencies in the optical variability properties reported 
herein and those found in previous analyses are proba- 
bly due to a biased sampling of high- redshift AGNs. Fur- 
thermore, compared to the DENIS samples, the UKIDSS 
samples reflect known properties of the optical variabil- 
ity. 

However, the near-infrared variability is negatively cor- 
related with the rest-frame wavelength, as has been sug- 
gested by previous studies. Nevertheless, it is notable 
that the near-infrared variability is anticorrelated with 
both the absolute magnitude and rest-frame time lag, 



which are opposite to the correlations with optical vari- 
ability. The anticorrelation with absolute magnitude can 
be explained by the near-inf rared emission from tori b y 
using the theoretical model of Kawaguchi & Moril (| 2 1 If) . 
The anticorrelation in the SFs should indicate that the 
near-infrared variability is more periodic over a longer 
time scale (up to ~10 yr). The detailed physical pro- 
cess behind the near-infrared variability needs to be re- 
vealed. By comparing the prope rties of our s a mples with 
the properties of sample s from lEnva et"aLl (|2002l ) and 
iNeueebauer fc Matthews! (|1999f ). we find that the gen- 
eral properties of the near-infrared ensemble variability 
are similar to the properties of radio-loud quasars. 

In addition, we examined the dependence of the near- 
infrared variability amplitude on optical variability, ra- 
dio intensity, and radio-to-optical flux ratio. Although a 
subset of our samples indicates that the optical variabil- 
ity is positively (negatively) correlated with variability 
amplitude at the J (K) band, these correlations might 
not be significant. The near- infrared variability ampli- 
tude tends to increase with increasing radio intensity and 
radio-to-optical flux ratio, and this positive correlation 
with radio intensity tends to be stronger in the range of 
F 20 cm > Q.01-0.1 Jy. 
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